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Disordered Quantum Solids 

T. Giamarchi, E. Orignac 



ABSTRACT Due to the peculiar non-fermi liquid of one dimensional sys- 
tems, disorder has particularly strong effects. We show that such systems 
belong to the more general class of disordered quantum solids. We dis- 
cuss the physics of such disordered interacting systems and the methods 
that allows to treat them. In addition to, by now standard renormalization 
group methods, We explain how a simple variational approach allows to 
treat these problems even in case when the RG fails. We discuss various 
physical realizations of such disordered quantum solids both in one and 
higher dimensions (Wigner crystal, Bose glass). We investigate in details 
the interesting example of a disordered Mott insulator and argue that in- 
termediate disorder can lead to a novel phase, the Mott glass, intermediate 
between a Mott and and Anderson insulator. 



1 Introduction 

Disorder effects are omnipresent in condensed matter physics, where one 
has to struggle very hard to deal with clean systems. Quite remarkably 
fermionic system exhibit marked differences with disorder in classical sys- 
tems. Indeed the very existence of a Fermi energy Ep "reduces" the effects 
of disorder since the relevant parameter now become the relative strength 
of the disorder compared to the Fermi energy D/Ep or the mean free path 
compared to the Fermi length kpl. Of course nature would not remain as 
simple as that, and quantum effects lead in fact also to reinforcement of 
the disorder effects and turn in low dimension a free electron system into 
an insulator, as pointed out by Mott and Twosejl). 

We have now gained a very good understanding of the properties of such 
disordered free electron systems. To tackle them an arsenal of methods 
ranging from diagrams^ |j| , scaling theory B , replicas^, |(J , supersymmetry M 
have been developped. Life become much less simple when interactions 
among fermions are taken into account. Although it is very intuitive to 
think that when interactions are small the noninteracting problem is a good 
starting point, such an intuition turns out to be wrong for a number of rea- 
sons: (i) even if in the pure system interactions can be "removed" from the 
system by resorting to Fermi Liquid theory^ this is not the case when dis- 
order is present. Because disorder renders electrons slowly diffusive slowly 
rather than ballistic, they feel the interactions much more strongly with 
explosive results M. Effective interactions increase when looking at low en- 
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ergy properties and Fermi liquid theory breaks down. The non-interacting 
physics is thus only relevant for very low disorder, (ii) When the dimension 
is small or the interactions strong to start with (like in systems undergoing 
Mott transitions) it is of course impossible to start from the non- interacting 
limit and one has to solve the full problem. What puts us at a disadvantage 
here is that most of the techniques useful for the noninteracting case also 
fail as soon as interactions are included: (i) supersymmetric method, which 
rests by construction on the quadratic nature of the hamiltonian is useless, 
(ii) it is now difficult now to really determine classes of diagram to sum. 

Disorder can still be averaged over by using replicas but then one is 
left with a complicated (and untractable) theory. Renormalisation group 
attempts have been made with some success [0, |n} but also with the prob- 
lem that the coupling constants diverge, so that the low energy fixed point 
remains elusive. Many other attempts in treating this very complicated 
problem exist in the literature and it is impossible to list them all[jl2j. Note 
that numerical studies are also hampered when studying this problem both 
because of the difficulty in taking the interactions into account (which more 
or less imposes either Monte-Carlo or exact diagonalization) and then to 
perform the complicated disorder averaging with enough statistics or large 
system sizes to get reliable results. 

A very peculiar situation occurs when one considers one dimension. On 
one hand one expects the difficulty to be maximum here. The interactions 
lead to very strong effects and destroy any trace of Fermi liquid giving rise 
to what is known as a Lutinger liquid. The disorder is also extremely strong, 
giving rise for the noninteracting case to a system so localized that the lo- 
calisation length is simply the mean free path and a diffusive regime is ab- 
sent. On the other hand one is in d = 1 in a much better situation to tackle 
the problem since the interactions can be treated essentially exactly using 
for example techniques such as bosonizationjl^, [l4|, [l5|, |l6|, [lTj so one only 
needs good techniques to tackle the disorder. Quite interestingly the physics 
of interacting and disordered one dimensional systems is the one of disor- 
dered quantum solids. Higher dimensional examples of such systems are 
the Wigner crystal, Charge Density Waves and Bose glass. In these notes, 
we will explain the techniques allowing to treat such systems, ranging from 
simple physical arguments to a quite sophisticated variational approach. 
In order to remain pedagogical and keep the algegra simple we mostly dis- 
cuss the technicalities on the simplest example of spinless fermions. We 
briefly discuss the specific physical realizations and give references so that 
the reader can look in more details the physical properties of these specific 
systems. 

Before we embark with the physics, let us point out that these notes 
results from the synthesis of various lectures. Some arbitrary choice of ma- 
terial had to be made in order to keep some level of clarity. Even if we have 
made some effort to cover various interesting topics, these notes cannot 
pretend to be as exhaustive as a full review. We thus apologize in advance 
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to anybody whose pet problem (or paper) is not covered in these few pages. 

2 Disordered interacting Fermions 

2.1 Model 

Let us consider spinless fermions hopping on a lattice with a kinetic energy 
t and an interaction V 

H = -tJ2(cU+i + h.c.) (1.1) 

i 

+ ^Vi-j{nz ~ n){n 3 - n) 
i>j 

where rij (n) is the local (average) electron density, and the rest of the 
notation is standard (t = 1 is the unit of energy). For nearest neighbor 
interactions this is the well known t — V model. When using a Jordan- 
Wigner transformation [jl8| to express the fermions in term of spins this 



latter model maps to an XXZ spin chain. In addition to (1.1) we want to 
submit the fermions to a disorder. We concentrate here on site disorder. 
Again most results/methods will carry over for randomness in the hopping. 
The randomness is simply 

Hint = ^ ~ ^- 2 ) 

i 

where fXi is a random variable of zero mean. 



2.2 Pure system 

The low energy properties of the pure system ( |l.l[) ar e by now well under- 
stood. We will thus review the bosonization of (1.1) only very briefly to 
fix the notations and refer the reader to the various reviews on the subject 
^4, 15, [u| [l7|. To get the low energy physics, it is enough to focus on 
the excitations aroung the Fermi points. In the continuum limit i — > x this 
amounts to express the fermion field in term of slow varying (with respect 
to the lattice spacing) field of right (with momentum close to +kp) and 
left (with momentum close to —hp) movers 



In term of these fields ( |l.l| ) becomes 



(1.3) 



H 



-ihvf 



dx(tjj\d x ip + (x)-^d x ip-(x))- 



dxidx 2 V (xi-x 2 ) p(xi) p(x 2 ) 
(1.4) 
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where the density reads 

p{x) = + iplip- + {e- l2kFX ^ + ip- + h.c.) 



(1.5) 



The remarkable feature in d = 1 is that all the excitations of the system 
can be reexpressed in term of the fluctuations of density. If one intr oduces 
a field cj>(x) describing the long wavelength part of the density (1.5) reads 



P{%) = Po 



tt (lira) 



(e 



-i2k F x+2c/)(x) 



h.c. 



where a is a lattice spacing, po = n/a and: 

± 

V0(r) = ^±4(x)<A±(z 



(1.6) 

(1.7) 
(1.8) 



9 is a similar field but associated with the long wavelength part of the 
current. (f> and Id = are canonically conjugate. Both the Hamiltonian 
and the fermion operator can be expressed in terms of these two fields 



H 

ip ± (r) 



h 1 ilr 



uK(V6(x)) 2 + -(V0(x)) 2 



,-i{±<f{r)-e{r)) 



(1.9) 
(1.10) 



Using ( |1.10[ ) in (1.5) one e asily recovers (1.6). For the free Hamiltonian (1.1) 
(with V — 0) one has in ( |l.9|) u = vf and K — 1. What makes the boson 
representation so usefu l is the fact t hat e ven in the presence of interactions 
the bosonized form of (LI) remains (|l.9| ) but with renormalized (Luttinger 
Liquid) parameters u and K. 



2.3 Disorder 



Similarly the disorder can be rewritten in the boson representation. It is 
natural to separate the Fourier components with wavevectors close to q ~ 
(rj) and q ~ ±2kp (£, £*). The disorder becomes thus 



Ha 



ti(x)[1f\.i/>+ + ipUj-} + (€(x)1>\.ilf- + h.c.) (1.11) 



where rj and ^ are now two independent random potentials, rj is real but 
when the system is incommensurate (2kp ^ tt) £ is complex. The potentials 
r\ and £ correspond respectively to the forw ard s cattering and the backward 
scattering of the fermions on impurities. (1.11) reads in the boson repre- 
sentation 



Ha 



-tj(x)—V^)(x) 

7T 



+ h.c. 



(1.12) 
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For incommensurate systems the forward scattering can easily be absorbed 
by a gauge transformation on the fermions. In the boson language this 
amounts to shift ^ </> — — J x dyri(y). Since £ is complex and hence has 
a random phase this shift does not affect the backscattering term (this will 
be even more explicit on the replicated Hamiltonian ( 1.46| ) below). Thus 
the forward scattering can be accounted for completely. Its main effect is 
to lead to an exponential decay of the disorder averaged density correlation 
functions. Since the current and the superconducting correlation functions 
depends only on 9 (or equivalently in an action representation on dt4>) they 
are not affected by the shift. Transport properties thus depends on the 
backscattering alone, as is obvious on physical basis since forward scatter- 
ing cannot change the current. For commensurate potentials the situation 
is mode complicated since £ is real (e +l2fcF = e~ t2kF ) and the forward 
scattering now does affect the backscattering term (in other words by com- 
bining a forward scattering and a scattering on the lattice one can generate 
a backscattering term) . This has consequences that we will examine in more 
details in section For the moment let us focus on the incommensurate 
case and get rid of the forward scattering. 



3 Tackling the disorder 



Disordered one dimensional particles are thus described by 



S/h 





1 


j dxdr 


2nK 



{d T (j>) 2 + v{dAf 



2nah 



h.c. 



(1.13) 



where we have rewritten (1.9) and (1.11) as an action. We have reintro- 
duced fi and other pesky constants to show explicitly the various physical 



limits. Note that although we arc mainly concerned here with fermions (1.9) 
describes in fact nearly every one dimensional disordered problem ranging 
from dirty bosons [|l9| 21 to disordered spin chains p2[, since all this 
problems have essentially the same boson representation. We will examine 
these different systems in more details in section [l| ( |1 . 13 ) emphasizes the 
physics of the problem. The electron system c an b e viewed as a charge 
density wavef23|, since the density varies as (|l.6|) 



p(x, t) = p cos(2kpx — 2<f>(x, r)) 



(1.14) 



The elastic term in ( 1.13 ) wants the phase of the density wave to be as 
constant as possible, and have a nice sinusoidally modulated density. The 
disorder term on the contrary wants to pin this charge density on the 
impurities by distorting the phase, as shown on Figure ^. The problem of 
localization of interacting fermions is thus very similar to the one of the 
pinning of classical charge density waves p5| . The charge density wave is 
here intrinsic to the one-dimensional interacting electron gas and not due 
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FIGURE 1. Impurities (black square) will pin the density wave. The phase <j> will 
have to change compared to the undistorted density (dotted line), which costs 
kinetic energy. 



to a coupling to phonons. The main features are nevertheless similar, the 
main difference being the fact that the effective mass of the "CDW" is 
much smaller in the absence of the electron-phonon coupling and hence the 
importance of the quantum fluctuations is much higher. Various quantities 
of special interest have a quite simple expression in terms of the bosons. 
The conductivity is simply given by 

i 

a(uj) = - [{d T <t>d T 4>) iu!n , q =Q] iUn ^ +iS (1-15) 

whereas the compressibility is given by a similar correlation function but 
with differents limits 

« = —o{dx4>dx<t>)iu„=o,q (1-16) 

The correlations of the 2Uf part of the density and of the superconducting 
order parameter ip+' l P- are respectively given by 

Xp ( r ) = ( e «*(r) e -«*(o)) (i.i7) 
Xs (r) = ( e «fl(r) e -i2 S (o)j (L18) 

where r — (x, r) and 7rn(a;) = \79(x). 

Given the highly nonlinear nature of the coupling to disorder, ( |1.13| ) is 
quite tough to solve. On the other hand, if £(x) was just a constant we 
would have a simple sine-Gordon theory for which a great deal is known. 
We will thus tackle the problem by increasingly sophisticated methods. 
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3.1 Chisel and Hammer 

In the absence of quantum fluctuations, (j) would be a classical field and 
we would have a good idea of what happens. This is the way Fukuyama 
and Lee p5[ looked at this problem. Such an approximation is of course 



very good for "phononic" charge density waves |23| since the quantum 
term is II 2 /(2M) and thus very small. For fermions this corresponds to the 
"classical" limit % — > 0, K — ► keeping K — K/Ti fixed, and thus to very 
repulsive interactions. In that case we can ignore all quantum fluctuations, 
and look for a static solution for <f>. It is of course crucial for the existence 
of such solution that the disorder does not depends on time. This solution 
4>q (x) describes the static distortion of the phase imposed by the random 
potential. In the absence of kinetic energy (V</>) 2 , it would be easy to 
"determine" <j)Q. If we write the random field £ as an amplitude |£(x)| and 
a random phase 2£, the disorder term writes 

dx\Z{x)\e l2{ ^ x) -^ x)) +c.c. (1.19) 

The optimum is thus for 4>q(x) to follow the random phase on each point. 
For point like impurities located on random positions Ri, |£| would just be 
the strength of each impurity potential and £ = kpRi- Thus (f>o(x) — C( x ) is 
the generalisation to any type of disorder (and in particular to the Gaussian 
disorder so dear to the theorist) of the physics expressed in Figure [l]: get the 
density minimum at each impurity. In presence of kinetic energy following 
the random phase would cost too much kinetic energy. We do not know 
exactly now how to determine the optimal 4>q (x) but we can do some scaling 
arguments. Let us assume that <f> remains constant for a lengthscale L\ oc . 
On this lengthscale </> takes the value that optimizes the disorder term, 
which now reads 



e ,:20 + c.c. (1.20) 



If we take for example a gaussian distribution for £ 



£(x)£*(x') = V5{x-x') (1.21) 

one gets because of the average of a complex random variable on a box of 
size Lioc that the disorder contributes as 



E dis = -^/VL^e 1 ^- 2 ^ (1.22) 

where S is some phase. It clear that the optimum energy is reached if <j>o 
adjusts to this (now unknown) phase. The global energy gain now scales as 
y/L\ oc . Between two segments of size L\ oc the phase has to distort to reach 
the next optimal value. The distortion being of the order of 27r the cost in 
kinetic energy reads 

£ki„ cx J- (1.23) 
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minimizing the total cost shows that the length over which <^q remains 
constant is given by 

£ioc oc (~) 3 (1.24) 

This tells us that the system does pin on the impurities and that below 
L\ oc the system looks very much like an undistorted system. Since at the 
scale L\ oc , cf> varies randomly the 2kp density density correlations will 
decay exponentially with a characteristic size L\ oc . It is thus very tempting 
to associate L\ oc with the Anderson localization length. Note that for the 



free fermion point L\ oc oc l/D instead of (1.24), so the above formula is 
clearly missing a piece of physics when K is not zero. Nevertheless from 
this simple scaling argument we have obtained: (i) the fact that classical 
CDW or very very repulsive fermions are pinned (localized) by disorder; (ii) 
the localization length; (iii) the fact that the ground state should contain 
a static distortion of the phase due to the disorder. Unfortunately we have 
no other information on </>o, which is certainly a drawback. 

Even with our limited knowledge of the statics we can nevertheless try 
to extract the dynamics. Let us assume that all deformations of the phase 
which are not contained in the static distortion are small and thus that we 
can write 

</>{x,t) = <t> (x) + 84>{x,t) (1.25) 

with 5cf)(x, t) <C 4>o (x) in a very vague sense since we deal with random 
variables. One can try to expand the random term in power of 5<fi 

S dis = J dTdx\£(x)\cos(2(<f>(x,T)-<;(x))) (1.26) 
~ -2 / dTdx\^x)\ cos(2(<Mz) - C(x)))(S^(x, r)) 2 



One can thus use in principle ( 1.26 ) to co mpute the various physical quan- 
tities. Note that the conductivity (|l,15j ) will not depend directly on the 
statics solution cf>o since dt4>o = 0, so we can hope to compute it. Of course 



the dependence of the fluctuations 5<f) in </>o is hidden in ( 1.26 ). If 4>q was 
following the random phase at every point, then the disorder term would 
just lead to a mass term for the fluctations and the optical conductivity 



would show a gap. In fact this is not true at every point so ( 1.26 ) leads to 
a distribution of masses for the fluctuations. Unfortunately the knowledge 
of cf>o is too crude to compute the conductivity accurately and depending 
on what exactly is L\ oc one can find either a gap, a non analytic behavior 
or a o~(ui) oj 2 behavior at small frequencies p5[ . Based on physical in- 
tuition Fukuyama and Lee opted for the later pM, but the method shows 
its limitations here and does not allow a reliable calculation of the physi- 
cal quantities. More precise calculations of (f>o and the conductivity can be 
performed in the classical limit K — > using a transfer matrix formalism 
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One thing that can be obtained from the partial knowledge of the fluc- 
tuations is the effect of quantum fluctuations pTj. Indeed when quantum 



fluctuations are present the expansion (1.26) is not valid any more since 
the cosine should be normal ordered before it can be expanded. This leads 
to 

cos(2(</> - C + $4>)) - -2cos(2(0 o - O)^ 2 ^ (W) 2 (1-27) 



where the average () has now to be computed sclf-consistently using ( 1.27 ) 
This leads to a modified localisation length p7fl of the form 



/ l \ 3-2K 

Lioc oc (- (1.28) 

This expression for the localization length suggests that a derealization 
transition is induced by the quantum fluctuations and occurs at K = 3/2. In 
the fermion language this corresponds to extremely attractive interactions. 



3.2 Starting from the metal: RG 

The previous method starts directly from the localized phase. It provides 
some limited information about this phase, but suffers from serious limita- 
tions. An alternative approach is to start from the pure Luttinger liquid and 
investigate the effects of disorder perturbatively, and build a renormaliza- 
tion group analysis. The RG provides us with the best possible description 
of the delocalized phase and the critical properties of the transition. It 
also gives a very accurate description of the localized phase up to length- 
scales of the order of the localization length L\ oc . Here again we describe 
the method for simplicity on spinless fermions and discuss more complex 
systems in section |[ 



To have a hint of the RG equations let us expand the disorder term ( 1.13 ) 
to second order. This leads to 

drdx f dT'dx'i(x)C{x')e 2l{<l}{r) ^ r ' )) (1.29) 



It is easy to see that at the tree level, (1.29) scales as T>L 3 2K . This leads 
to the scaling of the disorder 

dT> 

— = (3 - 2K)V (1.30) 
dl 

This traduces in fact the dressing of the scattering on the disorder by the 
interactions and has been derived using either diagrams or RG ps|, [ffl , 
|32| , |33| , |l9| , M . In itself it seems to confirm the result of ( [1.28] ) i.e. 
the existence of a transition. Note that the advantage of the bosonization 
derivation is to allow to reach the non perturbative point in interactions 
where such a metal-insulator transition would take place. 
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In fact (1.30) would not allow in itself to really determine the metal- 
insulator transition point. This can be seen by using the RG to compute the 
finite temperature (or finite frequency) conductivity of the system jl^. ^(J . 
The idea is simply to renormalize until the cutoff is of the order of the 
thermal length It ~ u/T corresponding to e l ~ It/ol. At this length scale 
the disorder can be treated in the Born approximation. As the conductivity 
is a physical quantity it is not changed under renormalization and we have: 



cr(n(0), D(0), 0) = a(n{l),D(l), I) 



n(l)D(0) e l D(0) 

g , n \nn\ = CT 0- 



n(0)D(Z) 



D(l) 



(1.31) 



where a(n(l), D(l), I) = a (I) and n(l) are respectively the conductivity and 
the electronic density at the scale I. a n — e 2 v 2 F /2TtKD is the conductivity 
in the Born approximation, expressed with the initial parameters. Using 
(1.30) one gets from (1.31) 



a(T) ~ —T 



2-2K 



(1.32) 



This result is the direct consequence of the renormalization of the scattering 



on impurities by interactions (1.30). One immediately sees that (1.32) alone 
would lead to a paradox since (1.30) gives a localized-delocalized boundary 
at K = 3/2 whereas (1.32) gives perfect conductivity above K = 1 (i.e. the 
nonintcracting point). One could also immediately see that if one introduces 
a new variable such as 

(1.33) 



V = e r al v 



the dimension of such a variable would be (3 — a — 2K) , leaving the location 
of the transition point as determined from ( 1.30| ) quite arbitrary. Although 
such a transformation seems arbitrary, if one considers that the disorder 
stems from impurities with a concentration rn and a strength V, the limit 
of Gaussian disorder corresponds simply in taking m — * oo with V — > 



keeping T> — riiV 2 fixed. Thus the choice a — 1 in (1.33) simply corresponds 
to T> = V 2 i.e. writting an RG equation for the impurity strength. 

The answer to this simple paradox is of course that fll.3CD should be 
complemented by another RG equation. In addition to renormalizing T> 
( 1.29 ) generates as well quadratic terms that renormalize the free part of 
the Hamiltonian, i.e. the velocity v and the Luttinger parameter K . Details 
can be found in |l^, 20 1. The main equation is the renormalization of the 
Luttinger parameter K and reads 



^ = -K 2 V/2 
dl 



(1.34) 



This equation describes the renormalization of the interactions by the disor- 
der. Both RG equations (1.3C) and (1.34) have a diagramatic representation 
shown on figure [| Using the flow ( 1.30] ) and (1.34) one can easily check 
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a) 1 a) 2 b) 

FIGURE 2. Diagrams describing the renormalization of the disorder by the 
interactions (a) and the renormalization of the interactions by the disorder (b). 
Solid and dotted lines are fermions with ±fej?, the wiggly line is the interaction 
and the cross is the impurity scattering. 




FIGURE 3. Phase diagram and flow for spinless fermions in presence of disorder, 
(a) is the flow in V and K variables, (b) the flow in the T> and interactions g. 
Disorder kills inelastic interactions. 
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that two phases exist as shown on figure |[ For large K one in the delocal- 
ized phase where the disorder is irrelevant and the system is a Luttinger 
liquid with renormalized coefficients u* and K*. All correlation functions 
decay as power laws, and because K* > 3/2 the system is dominated by 
superconducting fluctuations. On the transition line the exponent flows to 
the universal value K* —3/2 Below this line T> flows to large values, in- 
dicating that the disorder is relevant. This phase is the localized phase. 
This is obvious for physical reasons but can also be guessed from the exact 
solution known for the noninteracting line K = 1 (and any T>) which be- 
longs to this phase. As can be seen from ( 1.3Cl|) and ( |l.34| ) the transition is 
Berezinskii-Kosterlitz-Thouless (BKT) likcQ H in the K > ^ variables. 

In addition to the phase diagram itself a host of physical properties can 
be extracted from the RG. The simplest one is the localization length. 
One can use that for T>(1) ~ 1 the localization length is of the order of 
the (renormalized) lattice spacing ae l . The full determination needs an 
integration of both ( 1.30| ) and ( 1.34 ). Close to the transition the divergence 
of the localization length is BKT like (setting K = 3/2 + rf) 



Aoc 



(1.35) 



Deep in the localized regime, and for weak disorder, a good approximation 
is to neglect the renormalization of K in (1.30). A trivial integration of 
(1.30) then gives back (1.28). This we see that the SCHA calculation corre- 
sponds in fact, both for the phase diagram and for the localization length 
to the limit of infinitesimal disorder. 



Out of the RG one can also extract, using (1.31) the behavior of the 
temperature or frequency dependence of the conductivity. In the localized 
phase this can only be used up to the energy scale corresponding to the 



localization length i.e. -E p i n = ksT p 



TlUJ r 



v/L\ oc . Below this 



lengthscale another method than the weak coupling RG should be used. 
We will come back to that point in section ||. Here again, although the 
full flow should be taken into account one can get an approximate formula 
by ignoring the renormalization of K, which leads to ( 1.32 ). For K < 3/2 
(including the noninteracting point) any small but finite disorder grows, 
renormalizing the exponents and ultimately leading to a decrease of the 
conductivity, even if one started initially from K > 1. A very crude way of 
taking into accou nt bo th equations ( 1.3C ) and ( 1.34 ) would be to say that 
one can still use (1.32) but with scale dependent exponents (see MM, BCj) 



a(T) ~ T 2 - 2K ^ 



(1.36) 



This renormalization of exponents and the faster decay of conductivity is 
in fact the signature of Anderson localization. The equivalent frequency 
dependence is 

a(to) (xlo 2K - 4 (1.37) 
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Here again this formula break down below the scale o-> P i n ~ v / L\ oc which is 
the pinning frequency. Similarly below L\ oc correlation functions can again 
be computed using the RG, but of course the asymptotic behavior cannot 
be obtained. 



4 Other systems and RG 



Despite its limitations to physics above E p i n the RG is an extremely efficient 
method given its simplicity. It allows in addition a perfect description of 
the delocalized phase and of the critical behavior, something unatainable 
through the methods of section 3.1 and allow s for i ntere sting extensions. 

First let us note that the equations ( 1.3C ) and ( |l.34| ) also describe the 
case of a single impurity ]j3(| |37| . Indeed in that case one can go ba ck to 
the definition V = riiV 2 and take the limit — > 0, i.e. T> — ► 0. (1.34) 
shows that in that case K cannot be renormalizcd since a single impurity 
cannot change the thermodynamic behavior. Only ( 1.30| ) remains, leading 
directly to temperature dependence of the form ( 1.32| ), and a localized- 
delocalized transition at K — 1. More details on such a relation between 
the two problems and on the remaining open question can be found in p8j . 



Quite remarkably the set (1.3C-1.34) seems wrong. Indeed K naively 
depends on the (inelastic) interations. Perturbatively, for the pure systems 
K = 1 — V/(2wvf)- If one start for K = 1, i.e. for the non interacting 
system, it would thus seem from ( 1.34| ) that the elastic scattering on the 
impurities can generate inelastic fermion-fermion interactions. The solution 
of this paradox is hidden in the precise way the RG procedure is build. In 
order to have t he ela stic nature of the scattering on impurities, the time 
integrations in ( 1.29 ) should be dome independently for r and t' . When 
one performs the RG one intr oduces a cutoff and imposes |r — r'| > a. 
Thus a part is left out of |h2§) which is 



T> I dx drdr' p(x,t)p(x,t') 

J \t—t'\ <a 

^ 2T>a J dx J dr p(x , r) p(x , r) 



(1.38) 



which is exactly an inelastic interaction term. Thus in fact K contains not 
only the original inelastic interactions V but also a small correction com- 
ing from the disorder itself. In order to determine the flow for V it is thus 
necessary to take this small correction into account [[l9| |2(J which gives the 
flow of Figure ^-b. One thus sees that the elastic case V = indeed remains 
elastic and also that for spinless fermions, the perturbative flow seems to 
indicate that the inelastic interactions are reduced by the disorder. This is 
compatible with the physical image that one would get at strong disorder: 
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fcrmions localize individually and since the overlap of wavefunctions is ex- 
ponentially small, so is the effect of interactions. One could thus naively 
expect that below L\ oc the effect of interactions are strong but disappear 
above L\ oc . As we will see in section ^| the variational approach confirms 
this image. For fermions with spins (to be discussed below) we expect only 
the interaction in the charge sector to vanish. Interactions in the spin sector 
do remain and lead to random exchange. 

Of course many more physical systems can be studied by this method. 
This is the case for spin chains, than can directly be mapped onto spinlcss 
fermions [fL8|. A spin chain under a random magnetic field along z |22j is 
directly the problem that we solved in the previous section, with K being 
the anisotropy (K = 1 for an XY chain and K = 1/2 for an Heisenberg 
one). However although for fermion problems it is unlikely that one is at 
a commensurate filling this is the natural situation for a spin chain (since 
the magnetization is zero in the absence of external field the filling of the 
equivalent spinless fermion system is n = 1/2). A spin chain with random 
exchange is thus like a commensurate fermionic system with random hop- 
ping. We will come back to this peculiar case in section [?]. Quite remarkably, 
in one dimension, bosons|lj| lead to physics similar to fermions. This is 
due to the fact that in one dimension statistics cannot be separated from 
interactions. Interacting bosons can thus be represented by a bosonization 
representation quite similar to the one of fermions |59j. The density can be 
written as 

p(r) =p - p d x 4> + e ml ~ 27Tp0X - 2 ^ (1.39) 

n 

very similar to the fermionic form, while the single particle operator is now 

^{r)~pl /2 e l0 (1.40) 
(note the difference with the fermi onic operator). The hamiltonian is still 



described by the quadractic form ( 1.13 ). Now K — oo for noninteractin; 



bosons and K = 1 for hard core ones. There is thus a transition |19|, |20| 
between a superfluid state (for K > 3/2) to a Bose glass state where the 
bosons are localized by the random potential at K = 3/2. 

In a similar way one can of course treat fermions with spins. Equations 
are more complicated since they involve the charge and spin sectors, and 
we will not discuss the full physics here but refer the reader to |^0| . Let us 
just draw attention to one interesting consequence of the renormalization 
equation of the disorder for the problem with spins, which reads 



dV 
~dl 



— = {3-K p -K (T -g ll _)V (1.41) 



where gi± is the backscattering interaction between opposite spins. For 
a Hubbard type interaction gi± = U. For spin isotropic systems K a = 
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1 + <7ij_/(27TOf) and gi± is marginal with a flow 

-Six (!- 42 ) 



rt7 

For more general spin couplings either g\ ± — > and K a — > K*, or gu_ is 
relevant and a spin gap opens. The physics of the localization transition 
depends thus on the sign of the interactions with special physical conse- 
quences. The transition point move from K p = 2 (for infinitesimal disorder) 
for repulsive interactions to K p = 3 for attractive ones. For a Hubbard type 
interaction for which 1/2 < K p < 2 the derealization point can never be 
reached and the system is localized regardless of the strength of interaction 
and disorder. Another consequence is that in the presence of spin degrees 
of freedom the divergence of the localization length at the transition is not 
BKT like any more. One could naively think that in the repulsive case 
the <7ij_ term could be omitted and that the renormalization of the disor- 
der could be written (3 — K p — K*). For spin isotropic interactions this 
misses an important part of the physics. Indeed, for a Hubbard interaction 



K p = 1 — 17 /(2ttVf). Substituting in (1.41) leads (for the initial steps of the 
flow) 

dt TTVp 

whereas the incorrect substitution at the fixed point would lead to (l+^^r)j 
leading to quite different physics. ( 1.43Q implies that for Hubbard type in- 



teractions repulsive interactions make the system less localized [|40| than 
for attractive interactions, i.e. that the > ■ Similar effects exist 

for the charge stiffness and the persistent currents, i.e. for a system with 
spins the persistent currents are in fact enhanced by repulsive interactions. 
This counter-intuitive statement can be explained physically: interactions 
have two effects: (i) they tend to reinforce, when attractive the supercon- 
ducting fluctuations in the system. This screens disorder and makes it less 
effective. This is the only effect occurring for spinless fermions. (ii) when 
spin degrees of freedom exists, repulsive interactions also tend to make the 
density more uniform by spreading the charge. This makes it more difficult 
to couple to disorder. These two effect compete. This has several conse- 
quences, in particular for mesoscopic systems. Of course, for fermions, true 
derealization can only be achieved with attractive interactions reaching at 
least the nearest neighbor. 

Many other systems have been treated by such method, there are one- 
dimensional systems with long-range 1/r interactions, that lead to a pinned 
Wigner crystal EL doped spin 1 chains [p2| , fermionic 0] and bosonic 
p4| ladders, spin 1 chains in a magnetic disorder J45|| , spin ladders pq . 
Since we want to focus here on the methods we refer the reader to the 
above references for a detailed discussion of the physics of such systems. 
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5 A zest of numerics 

Although we are mainly concerned about analytical method in these notes, 
let us mention some numerical results and methods that have been used 
in connection with the RG predictions. Although numerical studies have 
become very powerful in one dimension for pure systems the presence of 
disorder complicates matters. Three main methods have been used. 

Exact diagonalizations, have been used to study both the phase diagram 



and the charge stiffness of both spinless fermions 47, |48| (or equivalently 



XXZ spin chains) and fermions with spins Ji9, 50 with short range or long 
range interactions J5lj]. Using the finite size scaling of the spin stiffness 

p s = x 9 q^P = f (xf - ) ' wnere E(9) is the ground state energy of the dis- 
ordered XXZ chain with boundary conditions — e lip S^ , the localization 
length L\ oc can be obtained KtJ, and is in good agreement with the RG 



results of section 3.2. The behavior of the correlation length close to the 
transition point appeared consistent with the predicted BKT-like behavior. 
The results also suggested that a finite disorder was needed to disorder the 
ground state for K > 3/2. A similar study with systems sizes of up to 
L = 18 sites was also made for XXZ spin chains with a random exchange 
in [p2~| . Analysis of persistent currents was also in agrement with the RG 
prediction of section ^. Unfortunately the exact Diagonalization approach 
of the last section is limited to zero temperature and small system size. 

In order to consider bigger system sizes, one can use Quantum Monte 
Carlo methods. In |53| , such a study was performed for disordered bosons. 
The superfluid density was obtained as a function of interaction for a given 
disorder strength. It was shown that for not too repulsive equations, there 
was a phase with a finite superfluid density. For more repulsive interactions, 
a phase with finite compressiblity by zero superfluid density was obtained, 
in agreement with the Bose Glass theory of section [l| 

The most promising recent method is the Density Matrix Renormal- 
ization Group. It been introduced in the recent years as a method spe- 
cially designed to calculate the ground state of correlated one-dimensional 
systems |54[] . This method has been also applied to the problem of the XXZ 
chain in a random magnetic field parallel to the z axis by Schmitteckert et 
al. j5j|. The authors of ]5j| have been able to consider system size of up to 
L = 60 sites, and average over several hundred realizations of the disorder. 
Localization and phase diagram were also in good agreeement with the RG 
predictions. 

Clearly, various numerical checks confirm the predictions of the RG. Un- 
fortunately so far only the phase diagram, stiffness and localization length 
have been computed. This is clearly related to the complexity of the prob- 
lem at hand. What would be extremely useful would be informations on 
quantities deep in the localized phase such as the single particle Green's 
function, the ac or dc conductivity. Analysis of such quantities would nicely 
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complement the RG analytical study and allow for comparison with other 
analytical techniques more suited for the localized phase such as the vari- 
ational method we analyze in the next section. 



6 Variational Method 



Let us now study this problem using a completely different and at first sight 
more formal method. As usual it is very convenient to get rid of the disorder 



from the start. Given the non quadratic nature of (1.13) supcrsymmetric 
methods are unapplicable and we have to turn to replicas. The idea of the 
replica method in itself is quite simple. If we want to compute an observable 
O we have to do both average over disorder and thermodynamic average 



(O) 



VVp(V)(0) v = / DVp(V) 



J D(f>e 



-Sv 



(1.44) 



The action is usually linear in disorder and for Gaussian disorder the dis- 
tribution of random potential is p(V) oc e / dx V ^ ; so the average would 
be quite trivial without the denominator in ( |l.44| ). The idea is thus to 
introduce n fields and to compute 



Dcjye 



(1.45) 



which is exactly the quantity we want to average over disorder in (1.44) 
if one takes the formal limit n — > 0. Since ( |l.45 ) has no denominator av- 
eraging over disorder is trivial. Of course there is a price to pay: before 
the averaging the replicas are all independent fields but the averaging in- 
troduces an interaction between them. We have thus traded a theory de- 
pending on a random variable V but a single field for a theory without 
disorder but with n coupled fields. Usually this is still a situation we are 
better equipped to solve because of the large number of field theoretic 
method dealing with "normal" (i.e. transitionally invariant actions). For 
the particular case 



1.13) the replicated action is 



S/h = f dxdr 



-{d T <t> a y + v{d x <t> a y 



(1.46) 



V 

(2na) 2 h 



/ dxdrdr' cos((f> a (x, t) - (f> b (x, r')) (1.47) 



where a = 1, . . . ,n is the replica index. Disorder averaging has coupled 
the replicas via the cosine term. Because the disorder is time independent 
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this coupling contains two fields that can be at arbitrary time and is thus 
highly non local. For fermions one usually prefers to go to frequency space, 
where this imply conservation of the frequency for each replica index, but 
this would not simplify things here because of the cosine. 

This is up to now a totally formal procedure and nothing has been ac- 



complished. (1.46) is totally equivalent to (1.13) and the difficulty is of 



course to solve it. Based on the RG equation (1.30) one could think naively 
that since the localized pha se cor responds to T> — ► oo it would be safe to 
expand the cosine term in ( 1.46 ). Unfortunately it is easy to check that 
fails seriously when n — > is taken (it of course works perfectly for a fi- 
nite number of field n > 2). In order to circumvent this problem let us try 
to improve over this simple minded expansion of the cosine. Let us try a 
variational ansatz. We introduce a trial action So 



s /h 



1 



2f3LK 



^2 ^2 <l>a(q, ^n)Gj{q,u n )(j) b (-q, -u n ) 



(1.48) 



where the propagators G 1 are our variational "parameters" . As usual 
rE^/T^r-Ifwe introduce 



(2tt) ' 



z 



V<t>e 



-s/h 



(1.49) 



We then have the variational theorem for the free energy F = —hlog(Z) 



F < F tI = F + (S- S ) So 



(1.50) 



Since Sq is quadratic, ( 1.50| ) can be in general computed quite explicitely 
as a function of the (unknown) propagators G. The "best" quadratic action 
Sq is thus the one that satisfies the saddle point equations 



dF tl 



dG ab (q,u> n ) 



= 



(1.51) 



which gives a set of integral equations allowing to determine the unknown 
functions G. 

The observables are simpl y defi ned by quantities diagonal in replica in- 
dices as can be seen from ( 1.451 ) . For some quantities such as the com- 
pressibility it is necessary to be more careful since one has to substract the 
average, which is usually zero in a pure system or after averaging over dis- 
order but non zero for a specific realization of the disorder. Let us introduce 
the various propagators (time ordering in r is always implied): 

B ab (x,T) = {[(j) a {x,T) -0 6 (O,O)] 2 ) = 

{G a a (0, 0) + G bb {0, 0) - 2G ab (x, r)) (1.52) 
G ab {q x ,uj n ) = (<f>a(qx,u n )cl) b (-q x> -u> n )) (1.53) 
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The compressibility is given by 

X {q,0J n ) = - / dx / dre-^ x -^ x 



x{T T {n(x,T) - (n(x,r)))(n(0,0) - (n(0,0)») (1.54) 

which leads to the average static compressibility \s = hm^o (lim^^o xili lu )) 
( see (p.. 16 )). When expressed in terms of the replicated bosonized operators 



(1.54) gives 



Xs = lim lim q 2 G c (q,ui) (1.55) 
where we introduced an important propagator: the connected one defined 

Without the replicas this method is nothing but the well known Self 
Consistent Harmonic Approximation (SCHA) , which is known to work very 
well for sine-Gordon type Hamiltonians. Such a method gives in particular 
correctly the two phases (massless and massive) . Extension of this method 
to disordered systems was done in the context of classical elastic systems 
such as interfaces |36|. In quantum problems another level of complexity 
occurs because of the aforementioned non locality of the interaction in time. 
But before going to these problems, specific to the quantum systems, let 
us illustrate the aspects of this variational method when applied to dis- 
ordered systems, on a technically simpler example (for which this method 
was extremely fruitful |57|): the case of classical periodic systems. 



6.1 A classical example 

Let us take the action ( |l.46[ ) but with only a single time integral for the 
disorder term. Such action would be the result of the average on a disorder 
both dependent on space and time. Of course such a disorder would be quite 



unrealistic for quantum problems. However (1.56) would be a perfectly 
natural Hamiltonian for a classical problem where z — vt is now just one 
of the spatial dimensions [ p8|, |57j . To make the analogy more transparent 
let us use z = vt, and replace the integral over x by an integral in d — 1 
dimensions. If denote by r the d-dimensional space variable r — (x, z) the 
starting action is 

S/h = J d d r^Y.^) 2 (1-56) 



V 



(2'Kct) 2 hv 

ab 



[ d d rcos(2(Mr)~Mr))) 



One can see that (1.56) is exactly the Hamiltonian describing an elastic 
system such as a vortex lattice or a classical CDW in the presence of point 
like defects in d dimensions, h plays the role of the temperature for the 
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classical system, the elastic constant c is given by c = 1/ (irK) with K = 
K/Tl and PqA/2 = Vhf (2-kol) 2 would be the correlator of the classical 
disorder @. 

If we call q — (q x ,uj) the ei-dimcnsional momentum, without loss of gen- 
erality, the matrix G~^(q) can be chosen of the form G^ 1 = cq 2 8 a b — o~ab- 
We obtain by minimization of the variational free energy the saddle point 
equations 

G- 1 (q)=cq\ a^^^e-^^ (1.57) 
(nay 



Using ( L5g -1.53) one obtains 

B ab {r) =hj -^L(G aa {q) + G bb {q) - 2 cos(qr)G ab (q)) (1.58) 

For this particular problem, the connected part is not affected by disorder. 
This is the consequence of a hidden symmetry (statistical tilt symmetry) 



of (1.56), whose disorder part is not affected by any local shift of a (r) 
such as (j} a (r) — * 4> a {f) + /( r )j where / is an arbitrary function. Such 
a symmetry does not exist for the time correlated disorder natural in a 
quantum problem, with important physical consequences on which we will 
come back. The only interesting equation here is thus the equation for 
the off-diagonal par t a a ^ b- Because of the locality of the interaction term 



between replicas in (1.56) the self energy a ab is simply a matrix of constants. 
Given the symmetry of the original action/Hamiltonian ( [1.56 ) by per- 



mutation of the replica indices, it is very natural to look for a variational 
ansatz with the same symmetry. This would mean that the G~ l matrix 
would have only (for each value of q) two independent values: the diagonal 
one G aa = G and the off diagonal one G a ^b- Such a matrix can easily be 
inverted for any n, and the analytic continuation for n — > gives 

G c = G-G a ^ b = -^ T (1.59) 
G c 

G a +b = ^ (1-60) 

^ (Gc ) 



Using these inversion formulas is it easy to solve for (1.57). In d > 2, 
B depends on G c only and thus o~ a ^b is simply a constant proportional to 
disorder. Given the gaussian nature of the trial action ( 1.48| ) the correlation 



functions such as the density-density can easily be computed 

WW)) = (Pa(r)pa(O)) cx e - 2S ^> (1.61) 
Using ( |l.5S| ) shows that G aa (q) ^V/q A leading to a growth 

B aa (r) ~ 2> 4 - d (1.62) 
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Although this solution is perfectly well behaved, a stability analysis of 
the replica symmetric saddle point shows that it is unstable. This can be 
checked from the eigenvalue A of the replicon mode (56l Bsl . 



*KD_ -mf^G c(p) f dS 2 
{ira) 2 J {2ir) d c 



(q) (1.63) 



A negative eigenvalue A indicates an instability of the replica symmetric 
solution. We introduce a small regularizing mass i n G c : G c (g) _1 = cq 2 + /j 2 
and take the limit jj, — ► 0. It is easy to see from ( 1.63j ) that for d < 2 the 



replica symmetric solution is always stable. In that case disorder is in fact 
irrelevant, due to the strong quantum (or thermal for the associated clas- 
sical system) fluctuations. For d = 2 the condition becomes /z 2 ^ -1 ) < 1 
for small /x. Thus there is a transition at K — 1 between a replica symmet- 
ric stable high temperature phase where disorder is irrelevant and a low 
temperature (glassy) phase where the symmetric saddle point is unstable. 
For the classical system this is the well known Cardy-Ostlund transition 
[[59) , with very interesting physical aspects of its own. This transition is the 
equivalent for time dependent disorder of the localization transition studied 



in section 3.2. For time dependent disorder one time integral drops in (1.29) 
and (1.30) would become (2 — 2K)T> giving the transition at K = 1. More 
details can be found in J58|, |6(], J57|. For 2 < d < 4 the replica symmetric 
solution is always unstable. 

One should thus look for another way of inverting the 0x0 matrices 
than the replica symmetric one. Fortunately such a scheme was invented 
in the context of spin glasses. Instead of having a single value a for the off 
diagonal term, one introduces a whole set of values. Let us briefly illustrate 
the procedure here, refering the reader to |6l| ] for details. Let us introduce a 
set of integers tuq — n, mi, m,k+i = 1 such that mj/mj+i is an integer. 
One cuts the matrix in blocks of size mi as illustrated on Figure |J Elements 
outside the blocks have the value uq. The procedure is then recursively 
applied for the inner blocks. At the last step the value on the block of 
size mfc is (Tfe and a diagonal value a. Quite generally the matrix is thus 
now parametrized by a diagonal element and a whole function (in the limit 
n — ► 0) cr(u), where u € [0, 1] (notice the range of variation of the nit when 
n — * 0). Such matrices can also be inverted in the limit n — > 0), albeit 
with more complicated inversion rules than for the RS solution. When one 
has a continuous function (see Figure ||) the RSB is said to be continuous. 
Simpler cases are of course a constant function (a single off-diagonal value) 
which is simply the RS solution or a function continuous by steps. We have 
represented in Figure || (see also Figure ||) the case of a two-step RSB. We 
denote G(q) = G aa (q), similarly B{x) — B aa (x), and parametrize G a b{q) by 
G(q,v) where < v < 1, and B a b(x) by B(x,v). Physically, v parametrises 
pairs of low lying states, in the hierarchy of states, v — corresponding to 
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FIGURE 4. Replica symmetry broken matrices (here for a 2-step RSB). Each 
pattern correspond to a different value <ji (see text), and the diagonal one (black 
line) is a. 
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states further apart. The saddle point equations become: 

2V 



= T^^ 2 ^ (1-64) 



where 

d d q 



B(0, v) = 2hj j^h(G{q) - G(q, v)) (1.65) 

B(0, v) corresponds physically to the mean squared phase fluctuations at 
the same point in space (r = 0) (for the associated classical system this 
bould be mean squared relative displacements of the same object) but in 
two replica states, or more physically in two different low lying metastable 
states. The large distance behaviour of disorder-averaged correlators is de- 
termined by the small v behaviour of B(0,v). We look for a solution such 
that cr(y) is constant for v > v c , v c itself being a variational parameter, and 
has an arbitrary functional form below v c . This corresponds to full RSB 
(see Figure |J). The algebraic rules for inversion of hierarchical matrices |OT 
give: 



d d q 2ha'{w) 
(27r) d (G c (g)-i + HH) s 



B(0, v) = B(0, v c ) + I dw I j^— - - \\, — 2 (1.66) 



where [cr](v) = ua(v) — J Q dw<j(w) and 



'Ifmw (L67) 



This is a simple numbe r. Ta king th e der ivative of ( 1.64 ) with respect to v, 
using [er]'(v) = va'(v), ( |l.66| ), and ( 1.64 ) again one finds 

Since the integral is ultraviolet convergent, we have taken the short-distance 
momentum cutoff to infinity, q is a simple number 

_ f d d g f 1 \ 2 _ (2-d)n 1 - d / 2 
° d ~ J (27r) d W + l) ~ 2 d + 1 sm{dir/2)r(d/2) ( ' 

with Cd=3 = l/(87r), Cd=2 — l/(47r). Dcrivating one more time one gets for 
the effective self energy: 

[a]{v) = {u/u Q f' e (1.70) 

where 6 = (d— 2) and vo = 8hcdC~ d / 2 / (4— d). The shape of [cr](u) is shown 
on Figure |^. The solution ( p..7C| ), is a priori valid up to a breakpoint u c , 
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(u) 



(b) 



FIGURE 5. Shape of the self energy a as a function of u. (a) is the full RSB 
solution occuring for 2 < d < 4. (b) a one step RSB (d = 2). 



above which [a] is constant, since <r'(u) = is also a solution of the varia- 
tional equations. u c can also be extracted from the saddle point equations 
and we refer the reader to J58| for details. The precise value of u c is unim- 
portant for our purpose but the existence of the two distinct regimes in 
[er](it) has a simple physical interpretation that we examine in section 6.2 
Using ( 1.70| ) one can now compute the correlation functions. Larger dis- 
tances correspond to less massive modes and is dominated by the small u 



behavior of (1.70). One obtains 



((0(r) - <X0)) 2 ) = 2h 



^(l"«»(«r))<3(ff) 



(1.71) 



G(q) = — (1 
cq z 



1 dv [a](v) Zd n „ 21 

v 2 cq 2 + [a}(v) ) ~ q d [ ' 



with Z d = (4- d)/(AhS d ) and 1/S d = 2 d - 1 ir d / 2 r[d/2}. Thus for 2 < d < 4 
this leads to a logarithmic growth, 



i 



-^d logbl 



(1.73) 



with Ad = 4 — d, instead of the power law growth ( 1.62 ) of the replica 
symmetric solution. Note that the amplitude is independent of disorder. 



6.2 If it ain't broken ... 

It is thus necessary to break the replica symmetry to get the correct asymp- 
totic physics (l/q d propag ator) instead of the 1/g 4 given by the replica 
symmetric solution. From ( |1.72 ) one easily sees that at large enough q (i.e. 
for short distances) one recovers the replica symmetric solution. Thus u c 
and [c(u c )] (see figure |^) define a lengthscale L ~ (c/Jfffuc)]) 1 / 2 above 
which the RS solution does not describe the physics. This lengthscale cor- 
responds to the pinning length similar to ( 1.24 ) obtained by balancing the 



1. Disordered Quantum Solids 25 



elastic energy with the disorder one. Here 

E e i ~ L d - 2 (1.74) 

E dis ~ (175) 

leading to the famous Larkin g ||§ length L ioc ~ (l/X>) 1 /(4-d) which 
corresponds to the distance for which relative displacements are of the 
order of the lattice spacing p3 (or for which the phase here is of the order 
2n) p5|. Below this lengthscale the system has a single equilibrium state. 



It can be described by simply expanding in the displacements in (1.13) 



H dis = J d d xf{x)cj>{x) (1.76) 

where / is a random force, derivative of the random potential V. It is easy 
to see that this model gives the 1/q 4 propagator. This breaks down when 
the displacements cannot be expanded, i.e. for distances larger than L\ oc . 
As anticipated by Larkin and Ovchinikov[|63| , pinning occurs that tends to 
keep a domain of size L\ oc in place despite the thermal fluctuations. 

This help us to understand the physics of the RSB solution. To illustrate 



it let us focus for simplicity to the case in d = 2. As can be seen from ( 1.70 ) 
by letting d — > 2 the RSB solution in that case is a one-step breaking as 
shown in Figure §. There are some additional complications but they are 
unimportant for our discussion here. A normal self consistent approxima- 
tion would approximate the energy by a simple Gaussian centered in <j> — 0. 
The only way to incorporate the pinning is thus to put a mass term in the 
propagator 

As can be immediately inferred from the discussion of section [Tl|, this is 
a much too crude approximation. The RSB solution is smarter and ap- 
proximates effectively the distribution of displacements by a hierarchical 
superposition of Gaussians centered at different randomly located points 
in space. A pictorial view of these two cases is shown on Figure |[ The dou- 
ble distribution over environment and thermal (or quantum) fluctuations is 
approximated as follows J5j . In each environment there are effective "pin- 
ning centers" corresponding to the low lying metastable states (preferred 
configurations). Since all q modes are in effect decoupled within this ap- 
proximation, for each q mode a preferred configuration (a state) is <fi a (q). 
They are distributed according to: 

P({Mq)})~Y[e~^ q2lMq)l2 (1-78) 



Each is endowed with a free energy f a distributed according to an expo- 
nential distribution P(f) ~ exp(w c //T) (here u c = T/T g ). Once these seed 
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(B) 

EM 




FIGURE 6. (A) The RS solution describes fluctuations around a single state, 
with a single spring constant. (B) the RSB solution takes correctly into account 
the fact that different regions of space correspond to different equilibrium values 
of <j>, with a priori different spring constants. The weight of a given configuration 
has been represented by the height of the shaded region. 
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states are constructed, the full thermal distribution of the q mode (f> q is 
obtained by letting it fluctuate thermally around one of the states: 

P^q) ~ J2 W a e-* ( « 2+L ro')l^- fc '«)l 2 (1.79) 

a 

where each state is weighted with probability W a — e~f"/ T / J2p eT^' T - 
One thus recovers qualitatively the picture of Larkin Ovchinikov as the 
solution of the problem with the replica variational method. The Larkin 
length naturally appears as setting the (internal) size of the elastically 
correlated domains. The full RSB case corresponds to more levels in this 
hierarchy of Larkin domains (in some sense there are clusters of domains 
of size larger than L\ oc ) and the way this hierarchy scales with distance 
reproduces the exponents for displacements and energy fluctuations. We 
refer the reader to |57) for more details on the use of the variational method 
in this context and for the physical properties of such systems. 



6.3 Quantum problems 

We now apply the same method to the quantum problem. Because of the 



non locality in time of the interaction in (1.46) the solution will have 
quite different properties. We again use the trial action (1.48), with the 
parametrization 

vG-£(q,u n ) = {{vq)2 ^ 8 ab - (T„6K) (1.80) 



Although, as in section p.l| a a b is still independent of q because of the 
locality in space it is now dependent on u> n . This would render the solution 
extremely complicated if it were not for a remarquable property of quantum 
systems. Off diagonal replica terms such as <J a ^b only exist for the mode 
uj n = The general argument |36], ^37j is that in each realization of the 
random potential V, the disorder does not depend on r. Therefore before 
averaging over disorder: 

G aby = (<Pa (x, r)(j> b (0,0)) = {<j> a (x, r)) (0,0)) = (0 a (x,O))(0 h (O,O)) 

(1.81) 

It is important to note that such a property crucially depends on the as- 
sumption that the hamiltonian is r-independent and on the fact that equi- 
librium has being attained. This is the case considered here. 

The static mode u> n — thus plays a special role. This is quite nat- 
ural in a time independent disorder. Quite naively one sees already that 
the properties of the variational solution will thus be very similar to the 
ones of point like (i.e. totally uncorrelated) disorder in d spatial dimen- 
sions (here d = 1). This suggest strongly that the variational method will 



pull out a static solution, reminiscent of the one introduced in section 3.1 
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and treat the fluctuations around this static solution. We will come back 
to this point later. The solution can be obtained quite in all dimensions 
and we refer the reader to p7] for details. We specialize here to the case 
d = 1. In this case two type of solution exist, a simple RS solution with 
(7^ = 0. This solution is stable for K > 3/2. It corresponds of course to 



the delocalized regime where the cosine term in (1.46) is irrelevant. The 



variational method correctly reproduces the (gaussian) delocalized regime, 
and the correct transition point, but of course misses the renormalization 
of the Luttinger parameters given by the RG. For K < 3/2 although an RS 
solution still exists it is unstable and physically obviously incorrect |37j. 
One should look for an RSB solution. In that case the correct solution is a 
one-step RSB solution (it would be full RSB for d > 2) of the type shown 
in Figure ||. 

dr(l — cos(u;„t)) 



G-^.O = — (vq 2 



2V 



— ) + 

v h(Tra) 2 



exp(-2hB(x = 0,r)) 



du exp(—2hB(u)) 



(1.82) 



with 



<r(q,uj n ,u) 



2Vv 
(7ra) 2 



f3exv{-?i2B{u))8 u . 



(1.83) 



These equations are still formidable to solve. A simple parametrization is 

1 

-2HB(t) 



vG c 1 (g,u„) 
2Vv ^ h 



:{{vq) 2 + col) + Ei(l - <y n , ) + I(u n ) (1.84) 



-2HB(u>u e ) 



(1 - cos(cj„r))dr (1.85) 



Si = u c (a(u > u c ) — <j(u < u c )) — [cr](u > u c ) (1.86) 

a(u) = ^L e -™ B Mp8 n , Q (1.87) 
(7ra) z 

The parameters Si, the breakpoint u c and the function I(u) n ) have to be 
determined self-consistently. Let us examine first the general properties of 
the solution. 



Since I(uj„ = 0) = Q it is easy to check from ( 1.16 ) that the compress- 
ibility is unchanged by the disorder, since the "mass" term Si goes also 
away at u n — 0. The variational method thus correctly reproduces that 
the compressibility of an Anderson insulator is still finite, and practically 
unchanged (for free electrons where we can compute it) from the value 
without disorder. Correlation functions are also easy to obtain. We just 
give here the important physical point without the explicit derivation |p^| . 
Because of the presence of Si in the propagator (1.84) they will be massive. 
This leads to 

B(x — ► oo , t = 0) — ► x/L ioc (1.88) 
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A full calculation of the correlations shows that B{t) grows until r ~ L\ oc 
when it saturates 

B(x = 0,r -> oo) -> Cste = (Po^x) 2 (1-89) 

Since the variational action is gaussian one has for the correlation of the 
2k f part of the density 

Xp (x,t) = e -2*(«,T) (1.90) 

Two different physical effects arc described by the above correlation func- 
tions. In the absence of disorder Xp( T ~ * 00 ) would go to zero as a power 



law, a sign of the wandering of the particles. ( 1.89 ) traduces the fact that 
the time independent potential localizes the particles at a given point in 
space, instead of letting them fluctuate (unboundedly) due to quantum 
fluctuations in the absence of disorder. This leads to density correlation 
in time going to a constant up to a Debye Waller like factor. The behav- 



ior ( 1.8SD thus shows that a static solution (f>o(x) exists. The variational 



approach thus provides a good justification for the static solution 0o on 
which the method of section || is built. The correlation length of this static 



solution is given by (1.88). Since this length controls through (1.90) the 
exponential decay of the spatial correlations of density it is thus related 
to the standard localization length. Note that here both L\ oc and the full 
static solution are determined by the variational method. Simple dimen- 



sional analysis on (1.84) shows that 

Lioc oc (1.91) 
The solution of the variational equations fl67f leads back to the expression 



( 1.28 ) for L\ oc . Quite interestingly ( 1.89 ) defines a length if one writes the 



constant in units of the fermion spacing. This length is the "width" of 
the fluctuations of the particles as shown in Figure 0. Note the difference 
between the two lengths l± and L\ oc . They are related through |67fl 

2 

ll = ^(tiK)ln(L loc /a) (1.92) 

At the transition one expects that L\ oc diverge s as L \ oc ~ cxp(b/ (K c — K) a ) 
with a = 1/2 from the RG (see equation ( 1.35| )). Thus relation ( 1.92 ) 
predicts that: 

l - ~ Jk^kW (L93) 

diverges as a power law, which could be measured in numerical simulations. 

Let us now look at the conductivity, for which we need I(oj n ). If we call 
I'(oj) and I"(u>) the real and imaginary parts of the analytic continuation 
of I(aj n ), then the conductivity is given by 

Ul" + Uj-U,* + I' +Sl) 

aM = (-a,2 + Sl+ J0 2 + (/») 2 ^ 
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FIGURE 7. Localization. 

It is easy to see that 3ta(uj) goes to zero at zero frequency and thus the 
phase is indeed localized. But because of the analytic continuation, the 
existence or not of a gap in the optical conductivity is not linked to the 
existence of a "mass" Ei but to whether I" is nonzero at small frequencies 
or not. The equation for I(u> n ) takes a particularly simple form in the limit 
K, h — > while keeping K fixed. In this limit, we can write 



where the scaling function / satisfies: 
fix) = 

It can be shown that for to - 



)■ 



(1.95) 



(1.96) 



In a similar way ( 1.94 ) gives 



(1.97) 



(1.98) 



Such behavior is in agreement with exact results in one dimension |g, [|] 
up to logarithmic prefactors. At high frequency, one can show that cr(u>) ~ 
uj 2K ~ a in agreement with the RG result. The resulting conductivity is plot- 
ted on Figure N. 
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CO 

FIGURE 8. The conductivity of the one dimensional Anderson insulator in the 
limit K -^0,h^0,K = f fixed. 



6.4 The fine prints 

The variational method is thus an extremely efficient method for this type 
of disordered problems. It allows to get most of the physical properties in 
the localized phase, and give the qualitative features of the transition and 
the delocalized phase. Its physics is very similar to the one described for 

The variational method determines the 
(but by taking into account the 



point like disorder in section |6.1 
"static" 



solution for the mode u> n 
effects of all modes) without being restricted to a single Gaussian. It can 
then correctly compute the fluctuations (both in q and w n ) around this 
solution, contrarily to the approximate method of section 3.1 which was 
impeded by the lack of knowledge of the static solution. 

In d = 1 an additional remarkable property can be seen. Although most 
of the properties are independent of the value of u c , the conductivity is 
strongly dependent on it. As for the static solution one has continuous 
replica symmetry breaking for d > 2 which goes continuously to the one 
step solution in d — 2. In these dimensions taking the value of u c out of 
the variational equations gives the correct conductivity (no gap in a(u>)). 
In d — 1, the u c obtained by minimizing the free energy would give an 
incorrect (gapped) conductivity. Another way to determine u c is to use the 
marginality condition ||S7f that corresponds to the instability of the replicon 
mode, similar to the condition (1.63). This condition, more dynamical in 
nature, coincides with the free energy value of u c for d > 2 but is different 
in d = 1. One can check that the marginality condition always gives the 
correct conductivity. Some arguments for why it is so were given in JoTj , 
but this point is not yet fully understood. This phenomenon has since been 
found to occur in other systems ]6q. 
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FIGURE 9. The gaussian approximation treats correctly the excitations where 
the phase stays in one of the minima of the cosine. They contribute to the a. c. 
conductivity at T — 0. Soliton excitations bring the phase from one minimum 
to another. With the cos(0) term these excitations only cost energy at the kink. 
When the potential is replaced by the quadratic approximation, the energy cost 
of the "wrong" minimim is very high, hence a poor result for the physical results. 
These excitations dominate the d.c. conductivity for T > 0. 

The success of the variational method is obviously linked to the fact that 
we are looking at small "quadratic" fluctuations around a certain (in our 
case highly disordered) solution. In this case the replacement 

_ COS (0)^I 2 (199) 

is very reasonable. Such an approximation is very good to compute corre- 
lation functions of the variable (f>. This is the case for the density-density 
correlation and the optical conductivity. What is missing in the approxi- 
mation ( |l.99| ) are the solitons that go from one minimum of the potential 
to another minimum. The energy cost of such excitations is grossly overes- 
timated by the variational method as is shown on Figure |^. 

This has important consequences for the calculation of correlation func- 
tions that involve the soliton creation operator e lS '. These correlation func- 
tions are found incorrectly to decay exponentially with distance at equal 
time and to be zero at unequal time. As a result, the variational method 
does not allow the calculation of Fermion Green's function (they involve 
the operators e^ 0± ^) nor superconducting correlations (which involve the 
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<|>(x) 



FIGURE 10. Hopping of charges that leads to a finite d.c. conductivity. This 
corresponds to soliton-type excitations in the phase variable <f>. 



operator e ). In the case of a disordered XXX spin chain, the situation is 
even worse, since S x oc cos(0), whereas S z ~ — + (—) x l a cos 20. There- 
fore, even in the presence of a disorder that preserves SU(2) symmetry 
(such as a random bond disorder), the variational method would lead to a 
spurious breaking of rotational symmetry. It would also give poor results for 
systems that include both the 9 and the field in the Hamiltonian. Exam- 
ples of such theories include disordered Hubbard ladders j43| , disordered 
spin ladders jl6| or XXX spin chain in a random fields |^2[ . 

The mishandling of soliton excitations also limits our knowledge of the 
transport properties at finite temperatures. Indeed the optical conductivity 
does not correspond to transport of charge but charge oscillations around 
the equilibrium positions, it is thus well described by our harmonic approx- 
imation as shown in Figure ^. On the contrary transport at finite T involves 
real charge displacements. Since p(x) ~ V</>, displacing a charge amounts 
to make a solitonic excitation in the field cf> as shown in Figure [k]. Indeed 
Mott's arguments, to compute the conductivity of noninteracting electrons 



in presence of phonons a(T) ~ exp 



- ) 



is strongly reminiscent of 
an instanton calculation as shown in Figure A similar argument applies 
to Efros and Schklovskii calculations J(3?|. An extension of the variational 
approach to finite temperatures would indeed still lead to a(T,u> = 0) = 
proof that it is missing the excitations that are important at finite tem- 
perature. Unfortunately no way to treat such solitons has been found at 
present despite some attempts |7(J. 

Despite these limitations, the variational method is up to now the only 
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analytical method giving information for such localized systems in the lo- 
calized phase. As with all variational approaches, some physical insight in 
the properties of the system under consideration is needed to determine 
whether the method as any chance of success. Clearly, one must apply this 
method only to systems that can be reasonably well understood qualita- 
tively from their classical action. Fortunately many systems fall in this 
category, and we examine some of those in the following. 



6.5 Higher dimension: electronic crystals and classical 
systems 

First the GVM can be used to study classical systems using the standard 
mapping t — ► z. The action ( |l.46| ) and its extension to higher spatial 
dimension describes elastic objects (lines in this case as shown in Figure p"i"| ) 
pinned by columnar (i.e. time or z independent) defects. This situation is 
realized for example in vortex in type II superconductors irradiated by 
heavy ions (creating the linear track of disorder). This system in 2 + 1 
dimension is equivalent to a d = 2 quantum bose system in presence of 
pins. In a similar way than in d = 1 (see section |6.3[ ) such system has 
a pinned phase (the Bose glass) |2l], [Hf. T he variational method can be 
used to describe the Bose glass phase p%. However contrarily to d = 1 
it cannot be used to go to the superfluid regime since to describe a two 
dimensional "melting" of the Bose glass phase dislocations are important 
(no dislocations exist in d = 1) and for reasons explained above the GVM 
overestimates the energy cost of topological excitations. Another way to say 



it, is that in d > 1 we loose the elastic description ( 1.39 ) of the Fermion 
or Boson operators. The GVM can thus only be used in phases where the 
particles are localized so that some elastic description can again be used. We 
can thus use the variational method in higher dimensions to study electronic 
crystals. This included Charge density wave, but also the two dimensional 
Wigner crystal of electrons. In such a phase the electrons are confined 
by their repulsion (and a in some systems a magnetic field). An elastic 
description can be used. Some level of quantumness is hidden in the elastic 
parameters ("size" of the particles, quantization of the phonon modes of 
the crystal). For such systems the calculation of the optical conductivity is 
particularly useful since it is one of the few probes of such systems. Since 
the physics of such systems would deserve a review of its own we will not 
dwell further on it here but refer the reader to 172], p73| for details. 



7 Commensurate systems 

When the filling of the fermion system is commensurate the physics dis- 
cussed above is modified in various ways since the backward scattering on 
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X 



FIGURE 11. The action ( 1.46 ) in d = 2 describes elastic lines in the presence of 
columnar pins. 



disorder becomes real. If the forward scattering still exists, not much is 
changed. Two special cases will thus occur: (i) the forward scattering is 
absent. This occurs because of a symmetry of the system. This is the case 
for example for fermions at half filling with a random hopping or for spin 
chains with random exchange, (ii) The commensurate potential (either due 
to the lattice or due to electron-electron interaction) would open a gap. 
There will be a competition between Mott physics wanting to get a com- 
mensurate (gapped) insulator and the disorder that would like to destroy 
such a gap (push the system locally away from commensurability) . 



7. 1 The peculiar random exchange 



For electrons at half filling with a random exchange the forward scattering 
does not exist and the disorder term is simply 



11 = I dxV(x)i[^)\.ip- -H.c.] 
V{x) ^ 



dx 



(27ra) 



sin(2^(a;)) 



(1.100) 
(1.101) 



Although this seems very similar to (1.13) one easily sees the difference 
on Figure 12. Contrarily to normal disorder where <j) follows the random 
phase of the random potential, here 4> = ±tt/4 depending on the sign of 
the potential. Thus 4> is nearly gapped but for its kinks. The low energy 
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FIGURE 12. The profile of <^(x) for a given realization of disorder V(x). As one 
can see, 4>(x) is only sensitive to the sign of V(a;) and not to its amplitude. As a 
result, <j> has kinks each time V(x) goes to zero. 

properties will thus be dominated by the kinks in </>. Such kink structure 
between the doubly degenerate minima makes it unlikely that the GVM 
can be used for this problem. 

In order to get an idea of the physics let use examine the case of free 
fermions K = 1. For E = one can easily construct the eigenstates by 
solving 



6V+ 
-dtp- 



- V{x)i)-{x) = 
+ iV(x)ip + (x) = 

to obtain the (unnormalizcd) 

ip+{x) ~ ifi_(x) > 



which obviously decays as 



ip±(x) 



(1.102) 
(1.103) 



(1.104) 



(1.105) 



If one wants to view such a state as an exponentially localized state this 
means that the localization length diverges when E — > 0. The divergence of 
the localization length is of course the signature that the physical properties 
will be quite different than the ones for the "normal" disorder studied in 
the previous sections. For the non-interacting case they have been studied 
using a variety of techniques such as replicas, Berezinskii methods, super- 
symmetry etc. [[74| [75|, f77|. The density of states diverge at the Fermi 
level asf7§ 

~ 1 n 1 ^ (1.106) 

Quite astonishingly if one considers that all states are still localized the d.c. 
conductivity is now finite [[74l [79| . This is surprising since it seems to violate 
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the Mott phenomena-logical derivation. The conductivity is proportional to 
the absorbed power. To make a transition between two localized states 
one need one occupied state, one empty one separated by an energy Too. 
The number of such states is p(e = 0)Tllu. The absorbed power (and the 
conductivity) is thus (up to log correction) 

o-(cj) oc p(e = 0){huj) 2 (1.107) 

correctly giving back for a constant density of states at the Fermi level 



uj dependence of section 5.3. One would thus naively have expected a 



(j{tjj) ~ u> in the case of a singular density of state (1.106). A qualitative 
understanding of the behavior of the d.c. conductivity is thus still lacking. 
When interactions are included, the problem becomes more difficult to 



tackle. Going back to the spin chain version of ( 1. 10C ) a real space renor- 
malization procedure has been introduced. This procedure works beauti- 
fully and allows the calculation of most correlation functions. We refer the 



reader to |8^, for details. As for the normal disorder (see section 3.2) 



interaction were found to be irrelevant at the disordered fixed point. The 
divergence of the localization length also changes drastically the correla- 
tion functions. In particular the spin spin correlation functions now decay 
as power law instead of being exponential. It would be extremely interest- 
ing to have an equivalent derivation of this real space RG directly in the 
boson representation. 

7.2 Mott versus Anderson 

A particularly interesting situation occurs when the non-disordered system 
possesses a gap. In that case the competition between this gap and the 
disorder is non trivial. This arises in a large number of systems such as 



disordered Mott insulators |82, [83|, |8j, |S5|, systems with external (Peierls 
or spin-Peierls system s p6|) or internal commensurate potential (ladders 
or spin ladders [|7| ||, |S9|,|9Cj], disordered spin 1 chains @ [42)). It is phys- 



ically simple to see that in order to destroy the gap one needs a disorder 
comparable to the gap (although in some specific cases Imry-Ma effects can 
destroy the gap for infinitesimal disorder |8^| ). This makes the complete 
description of the gap closure and of the physics of the resulting phases is 
extremely difficult with the usual analytic techniques such as the pertur- 



bative renormalization group described in section p.2| , due to the absence 
of a weak coupling fixed point. Indeed both the commensurate Mott insu- 
lator (which can be described by a sine-Gordon Hamiltonian [j92| ) and the 
disordered Anderson insulator correspond to strong coupling fixed points. 

To be concrete, but keep the technical details to a minimum let us con- 
sider again spinless fermions. The commensurability can be described by 



adding to the Hamiltonian (1.9) a periodic potential at the Fermi level 

dx(i/j\.ip-+ipt.il)+) (1.108) 
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which becomes in the boson representation 

H com = / dx-9--cos24> (1.109) 



The reason to consider here a periodic potential is specific to pathologies 
associated to the Mott gap for spinless fermions[^2j. Similar results are 
expected for a Mott insulator (which is due to the Akp potential of the 
lattice). 

If this system is studied by RG then one is faced with the competition 
of two strong coupling fixed points since 

§ = (2-K)g (1.110) 
^ = (3-2K)V (1.111) 

Using the usual qualitative argument consisting in taking the most diver- 
gent operator, it was concluded [^3| that if disorder reaches strong coupling 
(T>(1*) = 1 > <?(£*)) first, we will be in the Anderson phase, with a localiza- 
tion length Iq — ae l ~ (^) 3_2A ' , whereas if the commensurate potential 

reaches strong coupling first, we will be in the Mott phase with a corre- 

i 

lation length (or soliton size) d ~ . The phase transition between 

the two phases occurs for Iq ~ d. This picture relies on the important as- 
sumption that there is no other stable fixed point than the Mott insulator 
and the Anderson insulator. Even if it was so, it would not be possible to 
determine the phase boundaries, nor determine what type of critical point 
separate the Mott and the Anderson insulator. In order to make progress on 
these issues and obtain a more complete picture of dirty one-dimensional 
Mott insulators, one needs to solve the problem non-perturbatively f93||. 



The methods of section (6.3) are well adapted since the problem can be 



cast in the sine-Gordon form. 



7.3 Variational approach 

For the Mott versus Anderson problem, the variational action reads: 



E 



dxdr 
2nK 



V 



v{d x (j) a ) 



{d T 4>a) 



dxdr cos 2<f> a 



(2-Kah) 2 



^JdxJ^ drdr' cos (2(<f> a {x, t) - 4(i,r')))(1.112) 



The term cos20 a in (1.112) is responsible for the opening of a gap. We 
search for a saddle point with a form of the variational connected Green's 
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Function slightly generalized with respect to (1.84): 

vG- 1 (q t w n ) = -\={u 2 n + v 2 q 2 )+m 2 + Ei(l - <5„. ) + /(«;„) (1.113) 

7TA 



Where the parameter Si and the function I(uo n ) satisfy the equations (l.J 
and ( 1 . 8 5| ) . We give here the main steps of the solution and refer the reader 
t o p3 , 94fl for more details. The variational self-energy satisfies the equation 
(1.87). Finally, m satisfies the equation: 



Agv 



-2fiG(0,0) 



(1.114) 



The important physical quantities are simply given such as conductivity 
and compressibility are simply given by (1.15) and (1.55). 

Since we expect the physics to be continuous for small enough K (i.e. 
repulsive enough interactions), one can gain considerable insight by con- 
sidering (see sec. 3.1) the classical limit h — > 0, K — > keeping K = K/% 
fixed. In this limit one can solve analytically the saddle point equations 
([1.84fl-( |L87l) , ( |l.ll3| H |LTT4l ) and compute m, Si and I{w n ). The result- 
ing phase diagram is parameterized with two physical lengths (for K — > 0): 
The correlation length (or soliton size) of the pure gapped phase 



d = 



4ff3T - 1/2 



or 



(1.115) 



and the localization (or pinning) in the absence of commensurability length 



lo = 



(av) 2 

16X>7T 



1/3 



(1.116) 



Contrarily to the naive direct transition predicted by the extrapolation 
of the RG, we find three phases |93|] as shown in Figure [13[ Their main 
characteristics in term of conductivity and compressibility are summarized 
in Figure |l4|. 
Mott insulator 

At weak disorder we find a replica symmetric solution with Si = but 
with m ^ (m depends on the disorder), m is given by the simple equation: 



in 



■ exp 



VK 



1/2 



j,2 7r 3/2 v l/2 TO 3 



(1.117) 



It is convenient to work with the physical lengths d and Iq and the corre- 
lation length £ where: 



r = — — • 

(irKm 2 ) 



(1.118) 
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Mott 


Mott 
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Glass 
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d/l 

* Conductivity * Pseudogap * Pseudogap 

Gap i n conductivity in conductivity 



* Incompressible * incompressible * Finite Com P ressibilt y 

FIGURE 13. The phase diagram of a system with commensurability and disorder. 
d is the soliton size, lo the localization (or pinning length) 




FIGURE 14. Conductivity in the Mott Insulator (solid line) and Mott Glass 
phases (dashed line). Insert: density n versus the chemical potential /i. A is defined 
in the text. 
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One can then rewrite ( |1.117| ) as: 



1 

d 2 



£2 = TT (, M> 



1 

16 



(1.119) 



For l /d > \ (f ) 1 , ( gTTTg ) has a single physical solution. For lo/d < 

5 (t) 1 ^' ( EH1 ) nas n0 solution, which means that the Mott insulator 
becomes unstable. In fact due to another contraint (see below) the Mott 
Insulator becomes unstable at an even smaller disorder. 

Having ra/0 leads to zero compressibility k — 0. Disorder reduces the 
gap created by the commensurate potential and thus increases £ compared 
to the pure case. Let us now examine the equation for I(uj n ) giving the 
transport properties. An expansion around h = in equation ( 1.85 ) gives 
the self-consistent equation for I(u) n ): 



8Vv 
(ira) 2 



rtn 



G c (x = 0,t)(1 — cos(Lo n T))dT 



Introducing the scaling form (to be contrasted with (1.95)): 
I{ui n ) = m 2 f 



irKm 



(1.120) can be recast in the form 
f(x) = A 



1 - 



y/1 + X 2 + f{ X ) 



where: 



A 



AVK 1/2 v 



(1.120) 



(1.121) 



(1.122) 



(1.123) 



Let us note that for A = 2, the equation ( 1.122D reduces to the equation 
flL.96| ). For A > 2 ( p..l22| ) has no physical solution. Using f jl.l23|) , this 
condition becomes: 

j < r (1.124) 

For A < 2, there is a physical solution of (1.122) such that \mi x ^ r ± 00 f(x) = 
1 + A and for x < 1 , f(x) = 1 + ax 2 + o(x 2 ) with a = A/(2 - A). The 
conductivity of the Mott insulator can be obtained from / in the form: 



cr(w) 



(1.125) 



7T {l + f(tx)-X 2 ) 

where x = uj/u* and u* — u/£ is the characteristic frequency associated 
with the correlation length £. The conductivity cr(ui) is zero if: 



OJ 



l + A 



2/3 



(1.126) 
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A 0(03) 




CO 



FIGURE 15. The frequency dependant conductivity of the Mott Insulator 



The replica symmetric phase with to ^ has thus a conductivity gap 
TiLUc and can be assimilated to a Mott insulator (MI). The behavior of 
the conductivity is plotted on figure ([l5|). Close to the threshold, one has 
<t{uo) ~ (lu — LUc) 1 / 2 . However the gap in the conductivity Tilo c decreases 
when disorder increases, and closes for d/lo — 2e~ 1 / 4 ~ 1.557 (see Figure 
|l6|) . For d/lo — * 2e -1 / 4 , the conductivity gap vanishes linearly with yA 
For d/lo > 2e -1 / 4 the MI solution becomes unphysical even though the 
mass to remains finite at this transition point, i.e. the system remains 
incompressible. For stronger disorder one must break replica symmetry, 
as for the pure disorder case of section 6.3. Here, however two possibilities 
arise depending on whether the saddle point allows for to 7^ or not. In the 
presence of a breaking of replica symmetry, one extra equation is needed to 
determine the breakpoint. As in the case of the Anderson insulator, such 
equation i s pr ovided by the marginality of the replicon condition discussed 



in section |6.4| Two phases exists: 
The Anderson Glass 

For large disorder compared to the commensurate potential d/lo > 1.861, 
m = is the only saddle point solution. The saddle point equations then 



reduce to those ( 1.84 )— ( 1.87 ) . Thus, we recover the A nde rson glass with 



interactions of section (6.3). As we have seen in section (3.3), in such phase 
the conductivity starts as a(u>) ~ lu 2 showing no gap and the compressibil- 
ity is finite. 

In the Anderson glass phase, the disorder washes out completely the 
commensurate potential. The MI and the AG were the only two phases 
accessible by renormalization techniques |^3| . Within the replica variational 
formalism however, we find that an intermediate phase exists between them. 

The Mott Glass 
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FIGURE 16. The dependance of the conductivity gap with disorder strength in 
the Mott insulator 



For intermediate disorder 2e -1 / 4 < d/lo < 1.861 a phase with both Ei ^ 
and m ^ is obtained. We shall call this phase the Mott Glass (MG). 
We shall not discuss in full detail the one-step solution of the saddle point 
equation here, we will rather stress the salient features of our solution. First, 
as a result of the marginality of replicon mode condition, m 2 + Ei remains 
constant in the MI and MG as disorder strength is increased (9^, p4| . In the 



MG phase, I{to n ) is still of the form (1.121) but m is replaced by y/m 2 + Ei". 



The reduced self-energy f(x) satifies ( 1.122 | ) b ut with A = 2 in the whole 
Mott Glass phase. This implies that (see ( 1.96 )) that the a.c. conductivity 
of the Mott Glass is identical to the one of an Anderson glass. However, 
since m ^ in the MG, the system is incompressible (re = 0) like a Mott 
Insulator. Thus, the Mott Glass is a new glassy phase (since it has Replica 
Symmetry Breaking) with characteristics intermediate between those of an 
Anderson Insulator and those of a Mott Insulator. 



7.4 Physical discussion 

The existence of a phase with a compressibility gap but no conductivity 
gap is quite remarkable since by analogy with noninteracting electrons |]84f 
one is tempted to associate a zero compressibility to the absence of avail- 
able states at the Fermi level and hence to a gap in the conductivity as 
well. Our solution shows this is not the case, when interactions are turned 
on the excitations that consists in adding one particle (the important ones 
for the compressibility) become quite different from the particle hole exci- 
tations that dominate the conductivity. A similar situation is obtained in 
the case of the one dimensional Wigner crystal |)5| , which has the conduc- 
tivity of a perfect Id metal, o~(u>) oc 5(u>) but a zero compressiblity since 
X = lim g _^o This argument suggest that the difference in one-particle 
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(A) 



U 



V v 



(B) 



u 



V 



FIGURE 17. (A) Energy needed to add one particle. (B) Energy needed to 
make the particle hole excitations entering in the optical conductivity. Because 
of excitonic processes, when disorder is added the gap in the optical conductivity 
can close before the single particle gap. This phenonenon, leading to the MG 
phase, should occur regardless of the dimension. 



and two-particle properties is a consequence of the strong repulsion in the 
system. 

In addition to the variational method itself the Mott glass phase can 
also be obtained by two other independent methods. Higher dimensional 



extensions of the present problem, similar to the one made in section 6.1 



can be treated around four spatial dimensions using a d = 4 — e functional 
renormalization group method (totally different form the d = 2 RG). Such 



study confirms [B3, 94 the ex istence of the intermediate Mott glass phase. 



One can also analyse ( 1.112 ) for zero kinetic energy and obtain the MG 
phase |Q. Although we have done the derivation of the Mott phase for 
fermions in one dimension we expect its physics to survive into higher 
dimension. This can be seen by looking at the atomic limit (zero hoping) of 
an interacting fermionic system (in any dimension) . If the repulsion extends 
over at least one interparticle distance, leading to small values of K , particle 
hole excitations are lowered in energy by excitonic effects. For example for 
fermions with spins with both an onsite U and a nearest neighbor V the gap 
to add one particle is A = U/2. On the other hand the minimal particle- 
hole excitations would be to have the particle and hole on neighboring 
sites (excitons) and cost A p .h. = U — V, as shown in Figure [l7| When 
disorder is added the gaps decrease respectively |m| as A — > A — W and 
Ap.h. — > Ap.h. — 2W. Thus the conductivity gap closes, the compressibility 
remaining zero (for bounded disorder). According to this physical picture 
of the MG, the low frequency behavior of conductivity is dominated by 
excitons (involving neighboring sites). This is at variance from the AG 
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where the particle and the hole are created on distant sites. This may have 
consequences on the precise low frequency form of the conductivity such as 
logarithmic corrections. When hopping is restored, we expect the excitons 
to dissociate and the MG to disappear above a critical value K > K* . Since 
finite range is needed for the interactions, in all cases (fermions or bosons) 
K* < 1. In addition we expect K* < 1/2 for fermions with spins. One 
interesting question is the question of d.c. transport in the three phases, 
and whether the Mott glass has a d.c. transport closer to the Anderson 
or the Mott phase. Since the excitons are neutral, one simple guess would 
be that such excitations would not contribute to the d.c. transport. The 
d.c. conductivity in the Mott glass phase would thus be still exponentially 
activated just as in the Mott insulator. Of course more detailed studies 
would be needed to confirm this point. 



8 Conclusions 

Many disordered fermionic system can thus be successfully described by an 
elastic disordered theory. In one dimension, this situation is ubiquitous due 
to the importance of collective excitations. Most physical system, whether 
one starts with fermions, bosons or spins, can be represented in terms of 
bosonic excitations. In higher dimension such a description is valid in crys- 
talline phases such as a Wigner crystal in which the quantum particles arc 
strongly localized due to their interactions or charge density waves. Dis- 
order then leads to rich physical phenomena coming from the competition 
between the elasticity, wanting a well ordered structure and the disorder 
that distorts the structure to gain pinning energy. This leads to the exis- 
tence of many metastable configurations and to glassy properties. Using the 
various methods described in these notes, we now have a good description 
of the low energy excitations of such structures. This gives access to a host 
of physical properties such as the a.c. transport properties. 

Clearly one of the most important open questions is the the issue of 
topological defects in such structures. Indeed, such defects are needed to 
describe the melting of these cristalline phases and will be necessary to go 
to more "liquid" phases in which the statistics (fermionic or bosonic) of the 
particles will play a much more crucial role. In addition d.c. transport is 
obviously dominated by such excitations. Unfortunately so far the methods 
able to tackle the properties in the localized phase such as the Gaussian 
Variational Method cannot handle such solitonic excitations, so radically 
new methods will need to be designed to handle them. 
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